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Abstract 

The Goldstone-Brueckner perturbation theory is extended to in- 
corporate in a simple way correlations associated with large amplitude 
collective motions in nuclei. The new energy expansion making use of 
non-orthogonal vacua still allows to remove the divergences originat- 
ing from the hard-core of the bare interaction. This is done through 
the definition of a new Brueckncr matrix summing generalized Brueck- 
ner ladders. At the lowest-order, this formalism motivates variational 
calculations beyond the mean-field such as the Generator Coordinate 
Method (GCM) and the Projected Mean-Field Method from a pertur- 
bative point of view for the first time. Going to higher orders amounts 
to incorporate diabatic effects in the GCM and to extend the projection 
technique from product states to well-defined correlated states. 
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1 Introduction 

The mean-field approximation relies on a particular choice of the approxi- 
mate trial-state of the system, namely a Slater-determinant | <I>o } [1] * de- 
scribing the system as N independent particles. The mean-field approxi- 
mation can also be viewed as the zero-order approximation of the actual 
ground-state energy in a perturbative expansion written in terms of the 
residual interaction. 

In the case of nuclear structure, the mean-field approximation fails from 
the outset because of the strong repulsive core of the bare nucleon-nucleon 
interaction which leads to divergences when limiting the perturbative expan- 
sion to any order in the interaction. However, Brueckner [2] showed that the 
energy expansion can be reordered as a function of hole-lines number in the 
graphs, which amounts to an expansion in the density of the system. Such 
a reordering takes care of two-body short-range correlations induced by the 
repulsive core of the nucleon-nucleon interaction. It leads to an expression of 
the energy in terms of a renormalized interaction G. Using it, one can study 
the lowest-order approximation and define a meaningful mean-field picture. 
One can then evaluate the improvements achieved by including higher order 
diagrams since each of them is well-behaved. 

In many-body systems, several ways to go beyond the mean-field approx- 
imation exist, depending upon the physical situation of interest [1]. With- 
out particularly relating it to any perturbative expansion, one can improve 
the approximate trial wave-function of the system in a variational picture 

*In its generalized form taking care of static pairing correlations, the mean-field ap- 
proximation makes use of a state being a product of independent quasi-particles instead 
of independent particles. 
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in order to include correlations in the ground-state and calculate excited 
states. The projected mean-field method [3, 4] and the Generator Coordi- 
nate Method (GCM) [5, 6, 7] are such examples. The correlations considered 
in these two methods are those associated with large amplitude collective 
motions which are known to be very important in low-lying nuclear struc- 
ture. This is done by using a trial state written as a superposition of several 
non-orthogonal product functions referring to different mean-fields: 

m = E/"i $ o) • a) 

a 

Once such a trial state is given, its mean energy: 



T. a ,pf k P *f k a {<\H\n) 

can be minimized with respect to variational parameters. The mean-field 
approximation is recovered if a single coefficient is non zero; its self- 
consistent version corresponding to the minimization of the energy with 
respect to individual wave-functions. 

A problem is that going beyond the mean-field through Eq. 1-3, it is 
no longer possible to connect the variational calculation to any existing 
diagrammatic picture. Indeed, pertubation theories developed up to now 
can only deal with small amplitude correlations, these correlations being 
short-ranged or long-ranged [1] . The treatment of large amplitude collective 
motions as defined by Eq. 1-3 makes use of non-orthogonal vacua connected 
through specific infinite sums of particle-hole excitations which can hardly 
be identified in usual expansions. To make such a link is a question of 
interest since the systematic character of a perturbative expansion is a very 
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powerful tool in order to classify the correlations included or forgotten in 
a given variational calculation. An appropriate perturbative scheme would 
also be useful to identify the effective interaction which aims at renormalizing 
two-body correlations in the context of large amplitude collective motions. 

The aim of the present paper is to extend the Brueckner-Goldstone the- 
ory in order to motivate the methods embodied by Eq. 1-3 from a pertur- 
bative point of view. The lowest-order approximation of the new expansion 
provides an energy having the form as given by Eq. 3 and thus describes large 
amplitude collective motions^. Higher-orders introduce more correlations 
and of different nature than in the lowest-order. Namely, they correspond 
to the introduction of non-adiabatic effects in the coupling of individual and 
collective degrees of freedom. Moreover, the perturbative scheme allows to 
eliminate divergences arising from the strong repulsive core of the interac- 
tion by summing Brueckner ladders in the context of mixed vacua. It thus 
provides a microscopically defined effective interaction to be used in Eq. 3. 

The present work is organized as follows: we first remind some important 

features of the usual Goldstone-Brueckner perturbation theory in section 2. 

Then we develop the new perturbative scheme in section 3 where we define 

and study a generalized Brueckner G matrix. In sections 4 and 5, the 

link between the lowest-order approximation and the configurations mixings 

used in variational calculations is illustrated for the GCM and the projected 

mean-field method. The conclusions are given in section 6. 

t- We concentrate on perturbation theories written in a Hilbert space with definite par- 
ticle number. In other words, only normal contractions (c ! c) are non-zero. Such theories 
aim at treating the particle-hole channel of the interaction. 
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2 Standard Goldstone-Brueckner Theory. 



2.1 Ground-State Expansion. 

Let us start with the Hamiltonian of the N interacting nucleons. It consists 
of a kinetic energy term and a two-body force V, which includes a realistic 
nucleon-nucleon interaction and the Coulomb interaction: 

H = t + V . (4) 

We define a one-body Hamiltonian potentially breaking symmetries 
of H. Its ground-state and excited states are denoted | <J>q ) and | re- 
spectively. We make the additional hypothesis that | <I>q } is non-degenerate. 
One specifies this mean-field through: 

h° = t + T a = ^e Qn «ta„ , (5) 

n 

K |<£?) = £? |<&?) , (6) 

where the {a} single-particle basis is chosen as diagonalizing fifi. For our 
purpose, we assume the index i to be discrete and the Sf to be ordered with 
increasing value. The exact Hamiltonian can be related to h,Q through: 



H = h% + h? , 

where 

h" = V — T a , (7) 

is characterized as the residual interaction. 

In what follows, we make use of the Gell-Mann-Low adiabatic theorem [8] 
which states that an eigenstate of H is obtained when evolving | <I>q ) adia- 
batically from t = — oo to t = 0. Actually, no rule exists saying that this 
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eigenstate is necessarily the ground-state of the full Hamiltonian. We will 
consider that it is so [9]. This theorem explicitly reads as: 



In | 0q ) , the label a reminds from which unperturbed state | <3?q ) the 
actual ground-state comes from. The adiabatic evolution operator U" in 
the a interaction representation [9] is: 



C/ e Q (Mo) = j h ° t/n U e (t,t ) e~ ih ^l h 

= ihgt/h ~ift H(e,r)dr/h -ih%t /h 



(9) 



where U e is the evolution operator in the Schrcedinger representation. The 
derivation of Eq. 8 relies on two fundamental properties of quantum me- 
chanics [9, 10]: 

(1) Making an eigenstate of H(to) to evolve adiabatically from t = to 
through the evolution operator of the system drives the system to an 
eigenstate of H(t\) at t = t\. 

(2) It is possible to find an adiabatic time-dependent Hamiltonian H(e, t) 
which reduces to the independent particle Hamiltonian h$ at to = -co 
and which equals the full Hamiltonian H of the system at t\ = 0. 

In the present case, the auxiliary time-dependent Hamiltonian switching 
the perturbation adiabatically is given by: 

H(e,t) = K + e" e l*l hi , (10) 
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where e is a small positive parameter. At very large times, both in the past 
and in the future, this Hamiltonian effectively reduces to which presents 
a known solution, whereas at t = 0, H(e, t) becomes the full Hamiltonian of 
the interacting system. The limit e — > is taken at the end and the results 
are independent of e. Note that this method does not rely on the application 
of the evolution operator corresponding to the actual Hamiltonian H on a 
given state "at to = — oo". 

The normalization of the actual ground-state | 6q } is (<3?q | 0q) = 1. It 
is chosen through the denominator of Eq. 8 which exactly cancels discon- 
nected vacuum/vacuum diagrams in the numerator [11]. This is crucial as 
both the numerator and the denominator diverge in the limit e — > 0. 

Performing explicitly the time integrations contained in U", Goldstone [11] 
found a linked expansion of | 0q } allowing for a diagrammatic representation 
and reading as: 

@ 0> = E ( ga ~ h . h lj I *0) Linked , (H) 

where Linked means that | $q ) does not occur as an intermediate state in 
the diagrams. 

2.2 Brueckner Ladders. 

Following Brueckner [2], each V involved in a linked diagram of Eq. 11 
is to be replaced by a reaction matrix G a summing particle-particle (p-p) 
ladders in the {a} single-particle basis. This removes the strong repulsive 
core of the bare nucleon-nucleon interaction which makes the expansion in 
terms of V irrelevant. The replacement of V by G a is possible because 
each time a V interaction occurs in the original expansion of |©q), one 
can find the infinite set of identical diagrams except that a succession of 
V interactions connecting particle states occurs instead of the original V 
and this, before any interaction takes place somewhere else in the graph. 
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The corresponding diagrammatic content is shown on Fig. 1. This matrix 
satisfies a self-consistent equation of the form: 



G a (W a ) = V + V W J_ h „ G a (W a ) , (12) 

with t 

Q a = E I a P a p') ( a P a p' I • (1 3 ) 

The Pauli operator Q a acts in the two-particle space by excluding oc- 
cupied states in | <&q ) as intermediate states. The W a dependence of G a 
denotes that the in-medium interaction of two particles depends on the en- 
ergy of the others during the interaction. The G a (W a ) matrix elements 
together with the precise definition of W a are given in appendix A.l. Iter- 
ating Eq. 12 shows that G a takes into account the interaction between two 
particles to all orders of the potential in presence of the other nucleons. 

Let us now replace each h% by Gf (W a ) = G a (W a ) - T a in Eq. 11. All 
diagrams with two successive G a interactions joined through particle states 
are excluded from the sum in order to avoid double counting; the sum runs 
over linked irreducible diagrams [11]: 



l @ 0) = E ( ga ~ h% g ?(^)) \n)Linked,Irred. ■ (14) 

All graphs are generated in this way. This is nothing but a reordering 
of diagrams. Each term of the perturbative expansion is now finite and 
"well-behaved" even if the original interaction contains a strong short-range 
repulsion as in the nuclear case [12]. The Brueckner scheme has to be seen 
as the minimum recipe to get rid of the hard core problem. More elaborate 
treatments of the many-body problem through systematic renormalization 

*The two-body states | a p oy ) are non-antisymmetrized here. 
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Figure 1: G a matrix elements summing p-p ladders with respect to the 
vacuum | In any ladder, all V interactions (except possibly the first and 
last) are drawn between two up-going lines simulating particle states outside 
the Fermi sea. The rest of the diagram is not allowed to have any interaction 
between the first and the last interaction in the ladder. Points appearing 
in the present figure are superfluous here, but will become essential in the 
following. They will be defined in Fig. 3. 
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of propagators and vertices exist, using for instance the self-consistent Green 
function approach [13, 14]. 



2.3 One-body Potential T a . 

We do not discuss in detail the question of the most appropriate one-body 
potentials T a to be chosen for p-p, p-h and h-h matrix elements [14, 15]. Let 
us just recall few points. First, one can define a fixed, purely phenomeno- 
logical one-body potential r a . Such a choice suffers from its non-systematic 
nature. Another choice is to define T a in analogy to the Hartree-Fock po- 
tential with the bare interaction V replaced by the G a matrix. This is the 
Brueckner-Hartree-Fock (BHF) approximation. It reads: 

{a t \T a \ aj )= (a iak \G a (0)\ aj a k ,)p« kaki , (15) 

for h-h, p-h and p-p matrix elements, where the one-body density matrix p a 
associated with | <J>q ) is defined in any single-particle basis {q} by: 

f - (c \F\Cr) ~ ^l^l*?) (16) 

Pc s c r - {C s \p \C r ) - ($a|$a) ' ^ 

and is diagonal in the {a k } basis in such a way that Pa k a k , = $kh $k'h, where 
h denotes hole states. 

In Eq. 15, G a is taken on the energy shell (W a = 0). Such a definition 
allows to cancel a large number of diagrams in the wave-function. This is 
schematically shown on Fig. 2. Diagrams 2a and 2b correspond to insertions 
on vertices. They involve a G a interaction with passive unexcited states and 
a T a p-h interaction. Diagrams 2c embody insertions on hole lines. These 
three groups of diagrams have been shown to cancel thanks to the generalized 
time ordering [16, 17] although some of the G a interactions are originally 
not on the energy shell (W a ^ 0). On the contrary, insertions on particle 
lines cannot be systematically cancelled by any definitions of (a p \ T a \ ay) as 
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Figure 2: Cancelled diagrams in | 0q) in relation with the Brueckner- 
Hartree-Fock definition of the one-body potential r a . Dashed lines mean 
that the remaining part of the summed graphs are identical. 
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depicted on Fig. 2d. The two parts of these diagrams have to be calculated 
cxplicitely when appearing. The choice of the p-p matrix elements of T a has 
been of considerable controversy but we will not go into more details here 
and refer to Ref. [18, 19, 20, 21, 22]. 

Finally, let us say that the BHF choice for r Q is not the most appropriate 
one from the quantitative point of view [15, 23, 24]. A definition including 
higher-order diagrams is known to take care of the rearrangement poten- 
tial [25] which significantly improves density distributions in heavy nuclei as 
well as the agreement between theoretical single-particle and separation en- 
ergies [24, 26]. Indeed, higher-order graphs are cancelled with such a choice 
in a way that the zero-order is a better approximation for observables such 
as the one-body density matrix determining one-body observables. 

2.4 Energy Expansion. 

The Gell-Mann-Low expression for the energy Eq = £q + A£ff of the inter- 
acting ground-state | 0q ) is obtained through the scalar product of H | 0q ) 
with the unperturbed bra ($q| and reads as: 

A£ ° -i^o {n \U? (0,-00)1*?) • (17) 
The Goldstone-Brueckner linked/irreducible expansion deduced from it 

is: 



A£ = E<*o|G?(Wa) ( — I - Gt{W a )j \n)Linked/Irred. ,(18) 

where the sum runs now over all connected graphs with no external lines. 
Finally, a meaningful mean-field approximation in nuclear-structure is de- 
duced by considering the lowest-order in terms of G a in Eq. 18. This is the 
Brueckner-Hartree-Fock approximation: 
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£ n=u = + G?(0) 

= (^|t + G a (0)|<I> a ) (19) 

— ta h a h Pa h a h + n ^ a h a h , a h a h X^) Pa h a h Pa h ,a h , > 
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where G^ mapQnQ(j (0) are antisymmetrized matrix elements taken on the en- 
ergy shell. One can effectively see that Eq. 19 formally has the form of a 
mean- field energy for the effective "Hamiltonian" t + G a (0). 

We will not go into more details concerning the standard Brueckner the- 
ory and properties of the Brueckner G matrix. For extensive presentations, 
we refer to Ref. [2, 11, 15, 27]. Let us just mention that the renormalization 
of short-range correlations associated with the exchange of to vector mesons 
is not the only significant origin of in-medium effects and corresponding 
density-dependence. For instance, the exchange of pions should be impor- 
tant. In fact, the Brueckner summation originally introduced as the minimal 
recipe to remove the hard-core problem provides in-medium effects, not only 
from the short-range repulsive part, but also from the long-range part, i.e 
from the tensor part of the bare force [28, 29]. This is due to the depen- 
dence of the tensor force contribution to the energy on the Pauli operator 
Q a . This suggests that a part of the in-medium effects generated by pions 
is taken into account through the ladder summation (the ladders cover the 
whole energy range from low to high excitations energies) . Evidently, such a 
perturbation theory based on nucleon-nucleon potentials (possibly including 
three-body ones) neglects the dynamical mesonic degrees of freedom. One 
should go back to richer treatments of the many-body problem including 
these degrees of freedom explicitely in order to evaluate quantitatively their 
full in- medium effects [30]. 
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3 Extended Perturbation Theory. 



The previous section reminded the Brueckner-Goldstone expansion referring 
to a given mean-field unperturbed state | 3>q). As said in the introduction, 
variational calculations relying on the mixing of non-orthogonal product 
states cannot be motivated by such a perturbative scheme. We thus develop 
a perturbation theory based on a superposition of non-orthogonal Slater- 
determinants. A natural idea would be to evolve a linear combination as 
given by Eq. 1 from t = — oo to t = 0. To do so, one has to find a one- 
body Hamiltonian h™ tx , whose ground-state is the mixing of non-orthogonal 
Slater-determinants, and define an auxiliary Hamiltonian of the form: 

H(e,t) = h^ ix + e- e|t| (H-h^ ix ) . (20) 

However, it is not possible to find h™ lx for a general configuration mixing 
such as the one defined by Eq. 1. To avoid this difficulty, we use a different 
approach based on the arbitrariness of the unperturbed Hamiltonian /iq with 
respect to which the ground-state of the system is developed in the standard 
Goldstone picture. This arbitrariness, which is usually used to optimize the 
unperturbed state through the "best" choice of T a , is used here to superpose 
the expressions of the actual ground-state obtained starting from several 
unperturbed Slater-determinants | 

3.1 Ground-State Expansion. 

Let us choose a set of mean-fields, corresponding residual interactions, asso- 
ciated non-degenerate ground-state and excited states {h^ , hf , | <J>o ), \ ®f }}■ 
For instance, {Iiq } may correspond to a set of identical deformed mean-fields 
with different intrinsic orientations. Their deformation corresponding to the 
minimum of the associated BHF energy, each of them has a minimal density 
of individual states around the Fermi energy. This makes any particle-hole 
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(p-h) excitation to cost some energy and each ground-state | $q } to be non- 
degenerate. A second case of interest consists in defining {h^} through 
constrained mean-field calculations along a collective deformation path. In 
this case, single-particle level crossings may occur at the Fermi level for some 
values of the constrained deformation a. This makes the corresponding | <£q ) 
degenerate. Thus, one has to pick up a discrete number of a values along 
the deformation path for which this does not occur. 

The rationale of the present derivation is to make each of these vacua 
evolve independently from t = — ootot = 0as processed in section 2. 
The normalization ( <£q | 0q) = 1 eliminating disconnected vacuum/vacuum 
diagrams is chosen for each a. In this way, the | 0q } simply differ through 
their norm. For each a, the standard Brueckner scheme is performed in 
order to get the ground-state wave-function under the form given by Eq. 14. 
Finally, we perform an arbitrary linear combination of the states | 6q): 

I @o) = ^ fa | ©0 ) = f a ( ga _ foa ^1 ) I ^0 ) Linked/ Irred. ■ 

(21) 

As each [ 6q) is the ground-state of H with the eigenvalue Eq, j 60) is 
also the eigenstate with the same eigenenergy. This is true whatever the 
coefficients f a are. Actually, only the normalization changes from Eq. 11 
to Eq. 21. This last equation constitutes the starting point of our extended 
perturbative scheme. The linear superposition in Eq. 21 seems to be redon- 
dant and trivial but it provides richer approximate states at any order of the 
expansion. In particular, the lowest order provides with a superposition of 
non-orthogonal product states as a well defined approximation of the actual 
ground-state of the interacting system. In the following, we will focus on 
the perturbative expansion of the corresponding eigenenergy. 
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3.2 Energy Expansion. 



In the standard formulation, the starting point the energy expansion is 
Eq. 17 which gives the diagrammatic development of the correction A£q 
to the unperturbed energy £ff. It is obtained through the scalar product of 
hf | 0q) with the unperturbed bra ($q |. In our extended scheme however, 
there is no well defined unperturbed state and unperturbed energy since 
each of the vacuum | <£q) evolves independently from t = — oo before the 
linear combination is performed. To avoid this ambiguity, we keep the sym- 
metry between the bra and the ket and directly express the total energy Eq 
through the mean value of H in the ground-state | Bo). Using Eq. 21, we 
get: 

<e |ff|e > = £o(e |e > 



= E /W^if^rb?) 0[t (22) 

a,/3,n a ,nf3 \ ^0 n / 

+ ^ ] _ fra^l J \®o) Linked/ Irred. , 



where the norm is: 



(G |6 )= E fpfa (*qI I G l —Q 7b\ ( 7^ ) l $ o) Linked/ Irred. 

a,P,n a ,n V £ ~ K J V^O " ^0 J 

(23) 

The Goldstone diagrams corresponding to Eq. 22 and Eq. 23 are irre- 
ducible in terms of the G 13 and G a matrices. The dependence upon the 
energy parameters Wp and W a is omitted for convenience. In opposition to 
the energy expansion in the standard Goldstone-Brueckner theory, Eq. 22 
will provide explicit correlation diagrams associated to the kinetic energy 
operator. 
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In the above equations, the diagrammatic expansion is written as a sum 
over (a, (3) of several types of vacuum/vacuum diagrams ($q| ... |<J>q) in- 
volving different starting and ending vacua. The sum runs over k 2 of these 
terms if k different vacua |<I>q) are considered. Considering a particular cou- 
ple (a, (3), let us give rules to calculate the corresponding matrix elements 
but also to define and draw the new associated diagrams. 

First, all operators having a subscript (3 in Eq. 22 and 23 are written in 
the corresponding single-particle basis which is associated with dashed 
lines in the diagrams. All operators having a subscript a are written in the 
corresponding {a} single-particle basis which is associated with solid lines 
in the diagrams. The Hamiltonian in the middle of Eq. 22 is chosen for 
symmetry reason to be written in the left {(3} and right {a} single-particle 
basis: 



where Vp m p pan0lq are antisymmetrized matrix elements. 

Second, one has to insert complete sets of particle-hole excitations be- 
tween these operators while keeping linked diagrams only. These are particle- 
hole excitations with respect to |$q) around operators with an a subscript, 
and with respect to |3>q) around operators with a [3 subscript. In order to 
describe associated graphs, several types of propagators have to be defined 
depending on their localization in the diagrams. Some of them together with 
basic definitions are given on Fig 3. For example, reading the formula 22 
from right to left, which is equivalent to read diagrams from bottom to top, 
one can be convinced that the nucleons will first propagate from |$q) to 
(<E>q | in {a} individual states. Propagators involving different vacua will ap- 
pear when going through H because of the complete set of states referring 
to |$q ) inserted on its right and the one referring to |$q) inserted on its left. 
After that, they will propagate from |<£q) to (<J>q| in {/3} individual states. 




(24) 
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Conventions 



= 1 single-particle annihilation operator 

= 1 single-particle creation operator 

v =1 empty hole state 

/s =1 occupied particle state 

Full line ) labelled a refers to the non-perturbated ket 

Dashed line v labelled |3 - refers to the non-perturbated bra 

• and ■ m ^ mean that the nucleons propagate from l^"/* to "x^jjl respectively 
Non-zero propagators are typically of the form (always drawn with two arrows) 

<OppX> s *\ <*>»:> 

' . p p — - = O • > I - — - 



Figure 3: Basic definitions for drawing graphs between different vacua ap- 
pearing in the generalized energy expansion 22. Examples of mixed propa- 
gators involved in these graphs are given. 
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Graphically, we choose to specify the vacua entering the contraction in a 
propagator by a square for |$q) and by a circle for |<5q) at the extremities 
of the line (see Fig. 3). Finally, apart from the presence of two different 
single-particle basis in the graphs, the use of the Generalized Wick The- 
orem [31] instead of the standard one [32] and thus the apparition of new 
propagators, other rules for the calculation of diagrams are identical to those 
defined in Ref. [11]. Because of the non-orthogonality between particle-hole 
excitations built on two different vacua \%) and non-zero diagrams 

are more numerous than usual. 

3.3 Generalized Brueckner Ladders. 

We have now to treat divergences brought about by the V interaction coming 
from the energy operator H in Eq. 22. Their treatment has been separated 
from the one of the divergences appearing in the time evolution of the wave- 
function, which have been taken care of through the Brueckner matrices G a , 
because it is of different origin. 

Each diagram of Eq. 22 involving a V interaction can be associated 
with diagrams differing only by the fact that a G a interaction connects the 
original two nucleons state to V through two {a} particle states and/or that 
a G 13 interaction after V diffuses the two nucleons from {/3} particle states 
into the two nucleons final state. This allows to replace V by a G^' a ^ matrix 
summing generalized particle-particle ladders where nucleons propagate in 
two types of particle states referring to two different non-orthogonal vacua. 
In this way, all diagrams originating from Eq. 22 are generated and each 
graph given in terms of G^'^ is well-behaved as will be shown in the next 
section. The diagram content of the new reaction matrix is illustrated on 
Fig. 4. 

Analytically, G^' a ^ is expressed in terms of G a and G 13 through: 
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G^ a \Wp,W a ) = V + G p {Wp) — - — -5 V 

(25) 

= G ( '(W f3 )V- 1 G a (W a ) , 

the last identity being obtained by grouping terms and using the self-consistent 
equation satisfied by G a (W a ) and G^{Wp). 

The G^' a ) matrix elements together with the definitions of Wp and W a 
are given in appendix A. 2. G^ ,a ^ also satisfies a self-consistent equation 
which is given in appendix A. 2. Then, the energy expansion can be re- 
written in terms of G^' a \ The Eq. 22 becomes: 



<e |tf|e > = £ $f a <*J| ( g(-^-s) 

a,f3,n a , n/3 \ b h 0j 



+ G (/3 ' Q) y ga _ ha G l) \^o)Linked/Ir 



where diagrams connecting G a to G^'°^ through two {a} particle states 
and/or G^' a ^ to C 3 through two {/?} particle states are excluded in order 
to avoid double-counting. 

The reaction matrix G^' a ^ is the effective interaction removing the hard 
core problem within the framework of our generalized expansion as shown 
below. A remarkable property of this effective interaction is its dependence 
with respect to the type of TV-body matrix element in which it is inserted in 



20 



. (M) 




G 

'^TTYTYYTj' 

i v i 



/\ 

/\ G 



(X 




G 



m 


V 


'■ 


li 




ii 


/\ 




A 


a 


A 


G 


A 




+ + 



'■ v w 
■ ■ 




Figure 4: Diagrammatic representation of the G^'°^ matrix summing gen- 
eralized particle-particle ladders. 
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Eq. 26. Thus, for k states entering the superposition defined by Eq. 21, the 
present scheme requires the calculation of k different G a matrices referring to 
each vacuum plus k 2 matrices referring to each type of iV-body matrix 

element. It may be very time-consuming for quantitative calculations. 

As the aim of a Brueckner matrix is to take care of short-range corre- 
lations which should not be very sensitive to the single-particle basis used, 
the matrices G a should be almost identical for all a. However, BHF calcu- 
lations have shown that this is not such a good approximation (see Ref. [33] 
for instance). We will not go further concerning necessary approximations 
to implement this scheme for quantitative calculations since we stay on a 
formal level in the present work. 

3.4 Pair Correlated Wave-Function and G^' a ^ Matrix Prop- 
erties. 

In standard Brueckner theory, it is useful to illustrate the G a matrix prop- 
erties through its action on an uncorrelated two-body wave-function §: 

[ a p a q ) = al a\ | } . (27) 

This defines a correlated two-body wave-function | ) [27] through the 
equation: 

V\^ q (W a )) = G a (W a )\a p a q ) , (28) 
which can be re-written using Eq. 12 as: 

\^(W a )) = \a p a q ) + w ®°_ - V\^ q {W a )) . (29) 
5 Except for the definition of Q a , the two-body states are antisymmetrized. 
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It can be shown [12] that the content of the Brueckner summation is such 
that the correlated pair wave- function ip° „(W a , r, R) vanishes rapidly as the 
separation distance r between the two nucleons decreases within the range 
of the repulsive core. The exclusion of the nucleons from the region of the 
core in the correlated wave-function makes finite two-body matrix elements 
of the form: 

(a r a s \G a (W a )\a p a q ) = {a r a s \V \^ q (W a )) . (30) 

This is the physical reason why the standard expansion formulae 14 
and 18 in terms of G a provides well-behaved diagrams. Let us now study 
the properties of the generalized G^' a ^ matrix entering the new energy ex- 
pansion. As we have symmetrized the role of the bra and the ket in the 
expansion, the quantity to look at is no longer some correlated two-body 
ket alone but the two-body matrix element ( /3 r j3 s \ G^ ,c ^ \ a p a q ) itself. 

Using Eq. 25 together with the definition 28 for the pair correlated wave- 
function, one can write 

{(i r p s \G^ a \W a ,W p )\a p a q ) = {f5 r p s \G? (W^V' 1 G a {W a )\a p a q ) 

(31) 

= (^ s (W p )\V\^ q (W a )) . 

This identity shows that G^' a ) is well defined to renormalize two-body 
short-range correlations stemming from the repulsive core of the interac- 
tion. It also shows how the energy expansion 18 provides well-behaved 
diagrams. Comparing Eq. 31 with Eq. 30, one sees that G^'°^ renormalizes 
more correlations than standard Brueckner matrices since it correlates both 
the two-body ket and the two-body bra. This originates from writting the 
ground-state energy as the expectation value of the hamiltonian in the ac- 
tual ground-state and resumming the ladders accordingly, instead of taking 
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the scalar product of H | ©o) with a vacuum and resumming the ladders as 
done in the standard Brueckner theory. 

Given the above rules and renormalized interactions G a , and G^ ,a \ 
Eq. 26 together with the norm provide a perturbative formula for the ground- 
state energy. It is worth mentioning that all the exponents (np , fia) defining 
the orders of the expansion for each TV-body matrix element in Eq 26 and 
Eq. 23 are all taken equal n a = np = n , V (a, (3). This means that, at 
a given order, identical diagrams with respect to each unperturbed state 
will be dropped both in the ket and in the bra^. The order at which these 
expansions are truncated must be the same for ( @o|#| ©o) and ( ©o| ©o}- 

3.5 Lowest-Order Approximation. 

Let us consider the lowest-order approximation. It corresponds to retain 
the n = term only in Eq. 23 and Eq. 26. The approximate energy of the 
interacting system becomes: 

which is identical to Eq. 3 except that the starting two-body interaction V 
has been replaced by the G^' a ^ matrix on the energy shell (W a = Wp = 0). 
Consequently, we have achieved the goal to provide general configurations 
mixing calculations with a well-defined perturbative equivalent in terms of 
an effective interaction removing the hard-core problem. The corresponding 
diagrammatic picture is shown on Fig. 5. This lowest-order already con- 
tains long-range correlations through the mixing of non-orthogonal vacua. 
This suggests that the presently developed perturbation theory should be 
appropriate for systems being soft with respect to some collective degree 

^ Once the hard-core problem is solved, the meaningful order in the energy expansion 26 
is defined as the number of explicit G" interactions associated with each term originating 
from the ket (or the bra). 
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Figure 5: Lowest-order approximation for the energy in terms of the G^' a ^ 
matrix. Square and dot symbols are defined in Fig. 3. 

of freedom since the mixing of several mean-fields is an efficient way to 
take the associated correlations into account. Within this framework, the 
lowest-order corresponds to an adiabatic Born-Oppenheimer approximation 
justified when the individual degrees of freedom relax much faster than the 
collective one. Going to higher orders in the expansion allows to include 
diabatic effects in the coupling of these two kinds of degrees of freedom. 

A mean-field approximation is recovered by taking a single f ao coefficient 
to be different from zero. In this case, Eq. 32 contains a diagonal matrix 
element only and G^' a ^ reduces to a single term G^ a °' a °\ In comparison 
with Eq. 19, the effective interaction appearing at the mean-field level in 
the present case is different from the standard one. As shown in section 3.4, 
Q(ao,ao) generates two-body matrix elements including more correlations 
than G a ° does. In some way, our scheme is not a simple extension of the 
standard Goldstone formulae since we do not recover it for f a = 5 aao f ao . 
Looking at this particular case, one could think that the expression of the 
energy 22 together with the G^ ao ' a °^ matrix lead to some double-counting 
with respect to the standard scheme. The Eq. 22 being exact, this is simply a 
different way of resumming two-body correlations into the interaction when 
expanding the whole energy instead of some correlation energy. It may 
appear as an unnecessary complication when considering a single Fermi-sea 
but is the simplest way of superposing several of them to get the energy as 
given by Eq. 32 at the lowest-order in G 13 and G a . This was our main goal. 
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The matrix G^ ,a ^ depends upon the matrix element (3>q| ••• | $0 ) m 
which it is inserted. Thus, the energy at the lowest-order cannot be factor- 
ize into a form as given by Eq. 2 which is built from a defined state | ^k) 
(see Eq. 1). Within the perturbative calculation renormalizing short-range 
correlations, | VP*,) should not be considered as the corresponding approx- 
imate wave-function of the actual ground-state of the system from which 
other observables could be calculated together with bare operators. Actu- 
ally, this is only by analogy with a true variational calculation making use 
of the starting Hamiltonian that one would use | Vl/fc) as the corresponding 
approximate state of the system. The calculation done in this way for ap- 
proximate values of other observables will be of the same quality as for the 
energy at a given order only if short-range correlations associated with the 
corresponding operator are small. This is precisely the situation which has 
been encountered in the GCM and the projected mean-field method when 
using phenomenological interactions depending on the mixed local density 
p(P' a \R) [7, 34, 35, 36]. Indeed, while this interaction depends on the matrix 
element (<&q | ... | <I>q), the state | ^k) has been used to calculate observables 
such as probability transitions. We will thus consider this approximation. 

The reason why no perturbatively well defined approximate state can be 

extracted from Eq. 32 is related to the fact that some diagrams are summed 

in G^' a ^ incoherently with respect to a well-defined order in the expansion 

of the state | ©o) • We have actually defined a perturbative expansion for 

the mean value of an observable (energy) without taking care of the state 

of the system itself. It is worth noting that defining an approximate energy 

through the mean value of H = t + V in any approximate form of | 0o) 

would lead to divergent diagrams. Indeed, V could not systematically enter 

a ladder in this case. That is why ladder resummations has to be performed 

in the energy before any approximation takes place" . 

"The situation is similar when dealing with approximate forms of | 0o) in E = £q + 
A£q = (<E>o|#|6o) in the standard Brueckner-Goldstone theory. Then, | $ ( ™ ) is not 
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The same kind of situation is encountered in usual shell-model calcula- 
tions making use of effective operators in a limited valence space. In fact, 
this expresses the fact that the nuclear mean-field or beyond mean-field en- 
ergy in calculations making use of phenomenological effective forces is more 
to be seen as a functional (of local densities in general) than as the mean 
value of a two-body Hamiltonian in a defined state. We anticipate the exten- 
sion of the Skyrme functional beyond the vicinity of a single mean-field [37] 
which will be proposed in a forthcoming publication. This extension rely- 
ing on Eq. 32, a formal link is kept with the bare interaction. Contrary to 
Kohn-Sham philosophy [38], this strategy is by nature systematic and does 
not aim at incorporating all the missing correlations at the considered order. 

The scheme developed up to now is valid whatever the coefficients f a 
are. While they are of no importance as long as no approximation is done, 
they will influence the results at all orders of the expansion. Let us now 
study two applications of the general scheme associated with two different 
choices of these coefficients. 

4 GCM. 

Once the set of product states is chosen along a collective deformation path, 
the perturbative energy can be minimized with respect to coefficients 
at any order of the expansion. At the lowest order, the minimization of 
Eo =0 [f% , fa] is equivalent to the eigenvalue problem: 
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(33) 

rigourously defined to be the approximate state associated with the BHF energy. 
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written in matrix form because of the dependence of the effective "Hamil- 
tonian" t + G (/3 ' q) (0, 0) on the matrix element ($^| ... | 

These are the Hill- Wheeler equations encountered in the GCM [5, 12] 
written in the non-orthogonal basis {| ^q)} in terms of the appropriate ef- 
fective interaction. This result motivates the GCM from the perturbative 
point of view for the first time. Let us emphasize that we are in no way refer- 
ring to the Ritz variational principle here since, as already mentioned in the 
previous section, the Eq. 32 is not obtained through the expectation value 
of the actual Hamiltonian into an well defined approximate state. In par- 
ticular, the corresponding energy can be lower than the actual ground-state 
energy. Finally, minimizing the energy at some higher orders in the expan- 
sion motivates the introduction of diabatic effects in the GCM [4, 39, 40]. 

Within the present interpretation, the diagonalization of the GCM ma- 
trix provides several sets of "eigenvectors" {fa}, k = 0,1, 2, .. corresponding 
to different approximations of the ground-state energy 32. Everything hav- 
ing been written from the outset for the ground-state, no excited state is 
obtained here. The corresponding | form a set of orthogonal states where 
each of them is not orthogonal to the actual ground-state | ©o) which it ap- 
proximates. Using a simple argument, it is however possible to interpret the 
I ^fc), k > 1, as good approximations of the excited states of the system. 
Indeed, while | ^o) has a maximum overlap with | 6o), | ^k) which is orthog- 
onal to | ^o) w ih have a large overlap with the eigenstate | @k) corresponding 
to the same ordering of the energies. 

5 Symmetry Restoration. 

The second case of interest consists in concentrating on the approximate 
state rather than on the approximate energy of the system. The actual 
ground-state | 6q ) defined through Eq. 11 obeys the symmetries of the sys- 
tem whatever those broken by the product state | <3?q ) from which it develops 
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at t = — oo. However, a problem arises when diagrams are dropped in the 
expansion, since then, symmetry breaking may appear in the wave-function. 
For instance, if the n = term only is retained, one recovers the symme- 
try breaking mean-field solution | <&q). Another way to stress this point is 
to remember that usual perturbative corrections with respect to | <£q ) will 
not restore any broken symmetry unless an infinite number of diagrams is 
summed** [1]. 

In variational calculations, it is well-known that such an infinity can be 
taken care of through a linear superposition of symmetry breaking states 
| <I>o ) as given by Eq. 1. This is the projected mean-field method. According 
to the symmetry, it corresponds to a particular choice of the unperturbed 
states | <3>q) and mixing coefficients f a . Let us now show how our extended 
perturbation theory motivates such a method by restoring good quantum 
numbers in the standard Goldstone-Brueckner expansion. We examplify this 
through the restoration of rotational invariance. 

We define | $>{]) as the non-degenerate ground-state of an axially sym- 
metric deformed one-body Hamiltonian A set of rotated states around 
the y axis orthogonal to the symmetry axis is obtained through : 

| ^) = exp (i ™ J y ) | <&°) = R(a) | , (34) 

with a varying from — n to n as an integer. The Hamiltonian being invariant 

under this tranformation, the rotating operator R(a) commutes with it. 

| <&{]) being /i[]'s non-degenerate ground-state with the energy £q, | <1?q) will 

be the one of the rotating Hamiltonian /ig = R{a) /iq R){a) with the same 

eigenenergy. Thus, the adiabatic evolution operator is written in the a 

**In the particular case of rotational invariance, one could start from a spherical mean- 
field and use Bloch-Horowitz perturbation theory for degenerate unperturbed states [41] 
if dealing with an open-shell nucleus. As we are interested in relating the projected mean- 
field method to a perturbative expansion, we omit this possibility and continue with 
deformed mean-field unperturbed states. 
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interaction representation under the form: 

U?(t,t ) = e lh « t/h U e (t,to) e~ ih ^l n 

= R(a) e ih W h R\a) e~ 1 ^ Ha ^W* R{a) e -<to/n jR f ( 

= R(a) f7°(Mo) R\a) 
because the a dependence of the auxiliary Hamiltonian is: 

H a (e,r) = R{a) H°(e,r) R\a) . 
Consequently, the ground-state 21 can be re- written as: 



(35) 



f — > — ' 

a 



(<^|[/f(0,-oc)|<^) 

_ U?(0, -oo)|*8> 



If f a satisfies [12]: 



e -o (*0|E/?(0, -oo)|*8; 



f /0 = 



2/ + 1 

2n 



sin(7ra/n) cIqq (na/n) , 



(36) 



(37) 



where d£f is the Wigner function for (I,M,K = 0) quantum numbers [1], 
the operator: 



(38) 
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is the projector on total spin / and projection on the z axis M = 0". Thus, 
the Eq. 36 becomes: 



I e ) = Pio I e°) , (39) 

the projection being trivial if no approximation is done on \Qq). For an 
even-even nucleus, the projection gives for / ^ and the identity for 
/ = 0. 

Finally, the extended perturbation theory allows to restore broken sym- 
metries in the standard Goldstone-Brueckner expansion by projecting the 
approximate state. This formulation motivates the standard projected mean- 
field method from a perturbative point of view for the first time and gen- 
eralizes it to any order of the expansion. Moreover, the method allows 
for a coherent summation of Brueckner ladders through the use of rotated 
Brueckner matrices R(a) G° W(a) for each component of the approximate 
wave- function. 

Let us now study the corresponding expansion of the energy. As Pjq 
commutes with H and satisfies Pf = Pjq, the actual ground-state energy 
defined by Eq. 22 becomes: 

E " - (esift. | e8> ' <40) 

which includes a simple sum over a. In order to get the form 26, the diver- 
gences associated with the V interaction coming from H have to be regu- 
larized. The generalized Brueckner matrix doing so reads as: 



^ The projector for M / or for a triaxially deformed product state | $§) is more 
complicated but the calculation can be extended without major difficulty. 
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+ VR(a)^-^G (W )R\a) 

(41) 

and depends on the matrix element (0q| ... R(a) | 6{j) in which it is in- 
serted. Thus, the projection presents a simplification with respect to the 
general case since a single standard G a=0 matrix and k G^°' a ^ matrices have 
to be calculated explicitely. Once two-body correlations are summed in the 
interaction, the energy at the lowest order in G a=0 is written: 

F n=0 = Z n a= -n f a ° (*°0\ [* + C(°.")(0,0)] g(a)|*jj) U9) 
" £S=-„/S°<*8l*(«)l*8> ' 

It takes the form of the projected mean- field energy on spin 1 = with 
the appropriate effective interaction G(°' a )(0, 0). Two remarks still have to 
be done. 

In agreement with the discussion given in section 3.5, the correlated 
wave-function : 

*o) = Poo I $8) , (43) 

is taken as the approximate ground-state of the system corresponding to the 
energy 43. 

Even if the projection of | 0{j) for / ^ is zero, this is not necessarily 
true for an approximation of this state. For the same reasons of orthogo- 
nality as in the GCM, it is natural to use Eq. 42 and 43 for / / as good 
approximations of the energy and the wave-function of the excitated states 
| @i) of the system. 
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6 Conclusions. 



We have derived an extended Brueckner-Goldstone perturbative expansion 
for the ground-state energy of an interacting system. This scheme allows 
to recover at the lowest-order the energy of a mixing of non-orthogonal 
product states in variational calculation. The expansion is valid for systems 
interacting through a strongly repulsive interaction at short distances as it is 
written in terms of a generalized Brueckner matrix renormalizing two-body 
short-range correlations in the appropriate way. 

The main achievement of the present work is thus to provide a link be- 
tween the generator coordinate or the projected mean-field methods making 
use of (phenomenological) effective interactions and a perturbative expan- 
sion of the actual ground-state energy of the system. In particular, our 
method allows to understand the usual projected mean-field method as the 
lowest-order approximation of a symmetry restored perturbation theory. 

In addition, this work formally provides a well-defined effective interac- 
tion removing the hard core problem to be used the GCM or the projected 
mean-field method. In particular, the summation of particle-particle ladders 
is adapted to the superposition of non-orthogonal vacua characterizing these 
two methods. The link between the presently defined effective interaction 
and phenomenological interactions used in configuration mixing calculations 
of finite nuclei such as the Gogny or the Skyrme forces is the aim of a forth- 
coming publication. A new prescription for their density-dependence beyond 
the mean-field relying on the presently developed perturbation theory will 
then be tested on quantitative configuration mixing calculations in finite 
nuclei. Let us say that the feasibility of perturbative calculations using this 
extended Brueckner-Goldstone theory to deal with large amplitude collective 
motions may be questionable at the present time. This is a reason why we 
have planned in a near future to use it mainly as a tool to get theoretically 
grounded guidelines for density-dependent phenomenological interactions to 
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be used in beyond mean- field calculations. 
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A G a Matrix Elements and Starting Energy. 

An explicit matrix element of G a as defined by Eq. 12 is: 



In this definition, W a corresponds to the total excitation energy of the 
intermediate state before the G a interaction. In Eq. 44, we took non- 
antisymmetrized matrix elements for V which brings about the necessity 
to antisymmetrize those of G a at the end of the calculation. The other pos- 
sibility is to start from the outset with antisymmetrised matrix elements of 



B G^ ,a ' Matrix Elements and Starting Energies. 

The reaction matrix G^' a ^ satisfies a self-consistent equation generalizing 
the one for the usual G a Brueckner matrix: 



(44) 




V. 



G^ a \W p ,W a ) 



V + G^ a \W p ,W a ) 



Q 



v 
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(45) 

+ V Q \ B GW> a \Wp,W a ) 



- V 



G^(Wp,W a ) V 



Then, an explicit matrix element is 



1 \ - 

+ 9 E 



^ ^ a €a v ~^~€a a €a r ^a s 



1 \ - 

+ 9 E 



2 ^ ^ rn +e f3n -e f3t -ep u -Wf 3 

e t > e Pu >e F 



VprrAM G%>;l ras (W^W a ) V ar 
4 ^ g {zpm+zpn-zpt-zpu-Wp) (e ap +e aq -e ar 

e Pt' e Pu> e F 
e a.r i e ots ' >€< fr 

In this definition, W a corresponds to the total excitation energy of the 
intermediate state before the G^ ,a ^ interaction while Wp corresponds to the 
total excitation energy of the intermediate state after the G^' a ^ interaction. 
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